Search for: All records

The 230 GHz lightcurves of Sagittarius A* (Sgr A*) predicted by general relativistic magnetohydrodynamics and general relativistic ray-tracing (GRRT) models by the Event Horizon Telescope Collaboration have higher variabilityM_ΔTcompared to observations. In this series of papers, we explore the origin of such large brightness variability. In this first paper, we performed large GRRT parameter surveys that span from the optically thin to the optically thick regimes, covering the ion-to-electron temperature ratio under strongly magnetized conditions,R_Low, from 1 to 60. We find that increasingR_Lowcan lead to either an increase or a reduction inM_ΔTdepending on the other model parameters, making it consistent with the observed variability of Sgr A* in some cases. Our analysis of GRRT image snapshots finds that the major contribution to the largeM_ΔTfor theR_Low= 1 models comes from the photon ring. However, secondary contributions from the accretion flow are also visible depending on the spin parameter. Our work demonstrates the importance of the electron temperature used for modeling radiatively inefficient accretion flows and places new constraints on the ion-to-electron temperature ratio. A more in-depth analysis for understanding the dependencies ofM_ΔTonR_Lowwill be performed in subsequent papers.

 

In 1977, Blandford and Znajek showed that the electromagnetic field surrounding a rotating black hole can harvest its spin energy and use it to power a collimated astrophysical jet, such as the one launched from the center of the elliptical galaxy M87. Today, interferometric observations with the Event Horizon Telescope (EHT) are delivering high-resolution, event-horizon-scale, polarimetric images of the supermassive black hole M87* at the jet launching point. These polarimetric images offer an unprecedented window into the electromagnetic field structure around a black hole. In this paper, we show that a simple polarimetric observable—the phase ∠β₂of the second azimuthal Fourier mode of the linear polarization in a near-horizon image—depends on the sign of the electromagnetic energy flux and therefore provides a direct probe of black hole energy extraction. In Boyer–Lindquist coordinates, the Poynting flux for axisymmetric electromagnetic fields is proportional to the productB^ϕB^r. The phase ∠β₂likewise depends on the ratioB^ϕ/B^r, thereby enabling an observer to determine the direction of electromagnetic energy flow in the near-horizon environment experimentally. Data from the 2017 EHT observations of M87* are consistent with electromagnetic energy outflow. Currently envisioned multifrequency observations of M87* will achieve higher dynamic range and angular resolution, and hence deliver measurements of ∠β₂closer to the event horizon as well as better constraints on Faraday rotation. Such observations will enable a definitive test for energy extraction from the black hole M87*.

 

Images of supermassive black hole accretion flows contain features of both curved spacetime and plasma structure. Inferring properties of the spacetime from images requires modeling the plasma properties, and vice versa. The Event Horizon Telescope Collaboration has imaged near-horizon millimeter emission from both Messier 87* (M87*) and Sagittarius A* (Sgr A*) with very long baseline interferometry (VLBI) and has found a preference for magnetically arrested disk (MAD) accretion in each case. MAD accretion enables spacetime measurements through future observations of the photon ring, the image feature composed of near-orbiting photons. The ordered fields and relatively weak Faraday rotation of MADs yield rotationally symmetric polarization when viewed at modest inclination. In this letter, we utilize this symmetry along with parallel transport symmetries to construct a gain-robust interferometric quantity that detects the transition between the weakly lensed accretion flow image and the strongly lensed photon ring. We predict a shift in polarimetric phases on long baselines and demonstrate that the photon rings in M87* and Sgr A* can be unambiguously detected with sensitive, long-baseline measurements. For M87*, we find that photon ring detection in snapshot observations requires ∼1 mJy sensitivity on >15 Gλbaselines at 230 GHz and above, which could be achieved with space-VLBI or higher-frequency ground-based VLBI. For Sgr A*, we find that interstellar scattering inhibits photon ring detectability at 230 GHz, but ∼10 mJy sensitivity on >12 Gλbaselines at 345 GHz is sufficient and is accessible from the ground. For both sources, these sensitivity requirements may be relaxed by repeated observations and averaging.

 

We introduce a new library of 535 194 model images of the supermassive black holes and Event Horizon Telescope (EHT) targets Sgr A* and M87*, computed by performing general relativistic radiative transfer calculations on general relativistic magnetohydrodynamics simulations. Then to infer underlying black hole and accretion flow parameters (spin, inclination, ion-to-electron temperature ratio, and magnetic field polarity), we train a random forest machine learning model on various hand-picked polarimetric observables computed from each image. Our random forest is capable of making meaningful predictions of spin, inclination, and the ion-to-electron temperature ratio, but has more difficulty inferring magnetic field polarity. To disentangle how physical parameters are encoded in different observables, we apply two different metrics to rank the importance of each observable at inferring each physical parameter. Details of the spatially resolved linear polarization morphology stand out as important discriminators between models. Bearing in mind the theoretical limitations and incompleteness of our image library, for the real M87* data, our machinery favours high-spin retrograde models with large ion-to-electron temperature ratios. Due to the time-variable nature of these targets, repeated polarimetric imaging will further improve model inference as the EHT and next-generation (EHT) continue to develop and monitor their targets.

 

The centers of our Galaxy and the nearby Messier 87 are known to contain supermassive black holes, which support accretion flows that radiate across the electromagnetic spectrum. Although the composition of the accreting gas is unknown, it is likely a mix of ionized hydrogen and helium. We use a simple analytic model and a suite of numerical general relativistic magnetohydrodynamic accretion simulations to study how polarimetric images and spectral energy distributions of the source are influenced by the hydrogen/helium content of the accreting matter. We aim to identify general trends rather than make quantitatively precise predictions, since it is not possible to fully explore the parameter space of accretion models. If the ion-to-electron temperature ratio is fixed, then increasing the helium fraction increases the gas temperature; to match the observational flux density constraints, the number density of electrons and magnetic field strengths must therefore decrease. In our numerical simulations, emission shifts from regions of low to high plasmaβ—both altering the morphology of the image and decreasing the variability of the light curve—especially in strongly magnetized models with emission close to the midplane. In polarized images, we find that the model gas composition influences the degree to which linear polarization is (de)scrambled and therefore affects estimates for the resolved linear polarization fraction. We also find that the spectra of helium-composition flows peak at higher frequencies and exhibit higher luminosities. We conclude that gas composition may play an important role in predictive models for black hole accretion.

 

We investigate general relativistic magnetohydrodynamic simulations to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterizing the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions.

 

The Event Horizon Telescope (EHT) recently released the first linearly polarized images of the accretion flow around the supermassive black hole Messier 87*, hereafter M87*. The spiraling polarization pattern found in the EHT images favored magnetically arrested disks as the explanation for the EHT image. With next-generation improvements to very long baseline interferometry on the horizon, understanding similar polarized features in the highly lensed structure known as the “photon ring,” where photons make multiple half orbits about the black hole before reaching the observer, will be critical to the analysis of future images. Recent work has indicated that this image region may be depolarized relative to more direct emission. We expand this observation by decomposing photon half orbits in the EHT library of simulated images of the M 87* accretion system and find that images of magnetically arrested disk simulations show a relative depolarization of the photon ring attributable to destructive interference of oppositely spiraling electric field vectors; this antisymmetry, which arises purely from strong gravitational lensing, can produce up to ∼50% depolarization in the photon ring region with respect to the direct image. In systems that are not magnetically arrested and with the exception of systems with high spin and ions and electrons of equal temperature, we find that highly lensed indirect subimages are almost completely depolarized, causing a modest depolarization of the photon ring region in the complete image. We predict that next-generation EHT observations of M 87* polarization should jointly constrain the black hole spin and the underlying emission and magnetic field geometry.

 

It is well recognized that population heterogeneity plays an important role in the spread of epidemics. While individual variations in social activity are often assumed to be persistent, that is, constant in time, here we discuss the consequences of dynamic heterogeneity. By integrating the stochastic dynamics of social activity into traditional epidemiological models, we demonstrate the emergence of a new long timescale governing the epidemic, in broad agreement with empirical data. Our stochastic social activity model captures multiple features of real-life epidemics such as COVID-19, including prolonged plateaus and multiple waves, which are transiently suppressed due to the dynamic nature of social activity. The existence of a long timescale due to the interplay between epidemic and social dynamics provides a unifying picture of how a fast-paced epidemic typically will transition to an endemic state. 

Abstract We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87. 

Gravitational lensing near a black hole is strong enough that light rays can circle the event horizon multiple times. Photons emitted in multiple directions at a single event, perhaps because of localized, impulsive heating of accreting plasma, take multiple paths to a distant observer. In the Kerr geometry, each path is associated with a distinct light travel time and a distinct arrival location in the image plane, producing black hole glimmer. This sequence of arrival times and locations uniquely encodes the mass and spin of the black hole and can be understood in terms of properties of bound photon orbits. We provide a geometrically motivated treatment of Kerr glimmer and evaluate it numerically for simple hot-spot models to show that glimmer can be measured in a finite-resolution observation. We discuss potential measurement methods and implications for tests of the Kerr hypothesis.